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ABSTRACT 

Gravito-acoustic modes in the Sun and other stars propagate in resonant cavities with a frequency below a given limit known as 
the cut-off frequency. At higher frequencies, waves are no longer trapped in the stellar interior and become traveller waves. In this 
article we study six pulsating solar-like stars at different evolutionary stages observed by the NASA Kepler mission. These high 
signal-to-noise targets show a peak structure that extends at very high frequencies and are good candidates for studying the transition 
region between the modes and the interference peaks or pseudo-modes. Following the same methodology successfully applied on 
Sun-as-a-star measurements, we uncover the existence of pseudo-modes in these stars with one or two dominant interference patterns 
depending on the evolutionary stage of the star. We also infer their cut-off frequency as the midpoint between the last eigenmode and 
the first peak of the interference patterns. By using ray theory we show that, while the period of one of the interference pattern is very 
close to half the large separation the other, one depends on the time phase of mixed waves, thus carrying additional information on 
the stellar structure and evolution. 
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1. Introduction 

Solar-like oscillation spectra are usually dominated by p-mode 
eigenfrequencies corresponding to waves trapped in the stel¬ 
lar interior with frequencies below a given cut-off frequency, 
Veut. However, the oscillation power spectrum of the resolved 
Sun shows regular peak structure that extends well above 
Wut (e.g. Jefferies et al. 1988[ Libbrecht 1988[ Duvall et al. 
|1991| l. This signal is interpreted as travelling waves whose 
interferences produce a well-defined pattern corresponding to 
the so-called pseudo-mode spectrum ( [Kumar et al.||1990 l. The 
amount and quality of the space data provided by the SoHO 
satellite ( Domingo et al.|ri995 i allowed us to measure these 
High-Frequency Peaks (HIPs) using Sun-as-a-star observations 
( jGarcfa et al. 1998|l fro m GOLF ( [Gabriel et al.|[T995| l and 
VIRGO (Frohlich et al. 1995 i instruments. The change in the 
frequency pattern between the acoustic and the pseudo-modes 
enabled the proper determination of the solar cut-off frequency 
( [Jimenez|2006| l. 

Theoretically, the cut-off frequency approximately scales as 
g ^Jp|Tef^, where Teff is the effective temperature, g is the gravity, 
and p the mean molecular weight, with all values measured at 
the surface. Hence, the observed Veut can be used to constrain 
the fundamental stellar parameters. However, solar observations 
showed that Veut changes, for example, with the solar magnetic 
activity cycle ( Jimenez et al.|20fT] l. Therefore, the accuracy of 
the scaling relation given above needs to be discussed not only in 
a theoretical context but also taking into account observational 
constraints. 

For stars other than the Sun, the detection of HIPs is chal¬ 
lenging because of the shorter length of available datasets and 
the lower signal-to-noise ratio (SNR). Nevertheless, in this study 


we report on the analysis of the high-frequency part of the spec¬ 
trum of six pulsating stars observed by Kepler ([Borucki et al. 


20101 for which we were able to characterize the HIP pattern and 


the cut-off frequency. The time series analysis and the prepara¬ 
tion of the spectra are detailed in Sect. 2. In Sect. 3 we describe 
how to estimate v^ut and compare it with our theoretical expecta¬ 
tions. In Sect. 4 we perform a detailed analysis of the HIPs and 
interpret the results as a function of the evolutionary stage of the 
stars. Finally we provide our conclusions in Sect. 5. 


2. Data Analysis 

In this study ultra high-precision photometry obtained by 
NASA’s Kepler mission has been used to study the high- 
frequency region of six stars with solar-like pulsations. The stel¬ 
lar identifiers from the Kepler Input Catalogue names (KIC, 
Brown et al. 201 Ijl are given in the first column of Table 1. 


Short-cadence time series (sampling rate of 58.85s, Gilliland 
et al. [2010[l up to quarter 17 have been corrected for instru¬ 


mental perturbations and properly stitched together using the 
Kepler Asteroseismic Data Analysis and Calibration Software 
(KADACS, [Garcfa et al.|2mT] l. 

For each star the long time series (~4 years) are divided 
into consecutive subseries and the average of all the power 
spectral density (AvPSD) is computed in order to reduce the 
high-frequency noise in the spectrum. As in the solar case (e.g. 
[Garcfa et al.[[199^ [Jimenez et al.[[2()05[ [2006| l subseries of 
~4 days (4 x 1440 points, 3.92 days) were used because they 
are a good compromise between frequency resolution (which 
improves with longer subseries) and the increase of the SNR 
(which improves with the number of averaged spectra). For the 
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same reasons (as also done in the solar case), we smooth - with a 
boxcar function - the AvPSD over 3 or 5 points depending on the 
evolutionary state of the star and the SNR of the spectra. From 
now on when we refer to the AvPSD we mean the smoothed 
AvPSD. An example of this spectrum is given in the top panel of 
Fig.[T]for KIC 11244118. Because of the short length of the sub¬ 
series, it is not necessary to interpolate the gaps in the data (see 
for more details [Garcia et al.||2014) , so we prefer the simplest 
possible analysis. However, we have verihed that the results re¬ 
main the same when using series that were interpolated with the 
inpainting algorithm (Pires et al.|2015|l. 


Table 1. KIC name of each star, number of initial 4-day sub¬ 
series, number of hnally averaged subseries after removing those 
with a low SNR, and initial and hnal Kepler short-cadence quar¬ 
ters used in the analysis. 


KIC 

#Series 

#Averaged series 

Qi-Qf 

3424541 

292 

270 

Q5-Q17 

7799349 

292 

268 

Q5-Q17 

7940546 

244 

184 

Q7-Q17 

9812850 

292 

276 

Q5-Q17 

11244118 

292 

265 

Q5-Q17 

11717120 

292 

276 

Q5-Q17 




500 1000 1500 2000 2500 

Frequency (p-Hz) 


Fig. 1. Top panel; Smoothed averaged power spectral density of 
4-day subseries of KIC 11244118. The red line is the severe 
smoothing used to normalize the spectrum (NAvPSD) shown in 
the lower panel (see text for details). 


To avoid subseries of low quality and to increase the SNR in 
the AvPSD we hrst remove all subseries with a duty cycle be¬ 
low 50%. Then, for each subseries of each star we compute the 
median of the flat noise at high frequencies, above 2 to 5 mHz 
depending on the frequency of maximum power of the star. From 
the statistical analysis of these medians, we reject those subseries 
in which the high-frequency noise is too high. We have veri¬ 
hed that the selection of different high-frequency ranges does 
not affect the number of series retained. A detailed study of the 
rejected series show that there is an increase in the noise level 
when Kepler has lost the hne pointing and the spacecraft is in 
‘Coarse’ pointing mode. This is particularly important in Q12 
and Q16. A new revision of the KADACS software (Mathur, 
Bloemen, Garcia in preparation) systematically removes those 
data points from the hnal corrected time series. In Table we 
summarize the number of 4-day subseries computed for each 
star, the number of subseries hnally retained for the calculation 
of the AvPSD, and the total number of short-cadence Kepler 
quarters available. 


Fig|^ shows the place of the observed stars in a seismic HR 
diagram. The black lines are the evolution sequences computed 


using the Aarhus Stellar Evolution Code (ASTEC Christensen- 


Dalsgaard 20081 in a range of masses from M - 0.9Mq to 


M = 1.6Mo in steps of M = O.IMq with a solar composition 
(Zq = 0.0246). The target stars cover evolutionary stages from 
late main sequence stars to early red giants. 



7000 6500 6000 5500 5000 

Te,, (K) 


Fig. 2. Seismic HR diagram showing the position of the six 
stars analysed. The position of the Sun is indicated by the O 
symbol. The evolutionary tracks are computed using ASTEC 
( Christensen-Dalsgaard|2008 1. 


3. The acoustic cut-off frequency 

3.1. Observations 

The observed stellar power spectrum has different patterns in the 
eigenmode and pseudo-mode regions. In particular, the mean 
frequency separation, (Av), between consecutive peaks is dif¬ 
ferent below and above the cut-off frequency. In the hrst case, 
it corresponds to half the mean large frequency spacing, while 
above Vcut it is the period of the interference pattern. In addi¬ 
tion, a phase shift between both patterns appears in the tran¬ 
sition region. The frequency where the transition between the 
two regimes is observed corresponds to the cut-off frequency. 
Hence, we ht all the peaks in the spectrum (actually doublets of 
odd, f = 1, 3, and even, ^ = 0, 2, pairs of modes due to the 
small frequency resolution of 2.89yL(Hz) from a few orders be- 
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fore Vmax up to the highest visible peak in the spectrum. To take 
proper account of the underlying background contribution (due 
to convective movements at different scales, faculae, and mag¬ 
netic/rotation signal), we prefer to divide the AvPSD by the same 
spectrum heavily smoothed (see the red line in the top panel of 
Fig. instead of using a theoretical model (e.g. |Mathur et al. 


|201 l| l. Although it has been demonstrated that a two-component 
model usually hts properly the background of stars (see for more 
details Kallinger et al.||20T4 l, the accuracy in the transition re¬ 
gion between the eigenmode bump and the high-frequency re¬ 
gion dominated by the HIPs is not properly described, which 
is the region we are interested in. Therefore, we use a simpler 
description of the background based on the observed spectrum 
itself. The length of this smoothing varies according to the evolu¬ 
tionary stage of the star. An example of the normalized resultant 
AvPSD (NAvPSD) is given in the bottom panel of Fig. [T] 

A global fit to all the visible peaks in the spectrum is not pos¬ 
sible for the following reasons; 1) it requires a huge amount of 
computation time owing to the high number of free parameters in 
the model, 2) the background is not completely flat around Vmax, 
and 3) the large amplitude dispersion of the peaks between Vmax 
and the HIP pattern biases the results of the smallest peaks by 
overestimating their amplitudes. We therefore divide the spec¬ 
trum into a few orders before Vmax and the photon-noise domi¬ 
nated region into two regions. We arbitrarily choose a frequency 
where the amplitudes start to be very small and treat separately 
the power spectrum before and after this frequency. We denote 
by p-mode region the low-frequency zone and pseudo-modes re¬ 
gion the other one. Several tests have been done in which the 
frequency separating these two regions were varied and in all 
cases the results remain the same. 

In the p-mode region we ht groups of several peaks at a time 
for which the underlying background can be considered flat. This 
ht depends on the evolutionary stage of the star and the SNR. 
We have checked that htting a different number of peaks in each 
group does not change the hnal result. Appendix A shows the 
analysis of all the stars, and in the caption of each hgure we ex¬ 
plicitly mention the number of peaks htted together. The pseudo¬ 
modes region can be htted at once because 1) the background is 
Hat (see the bottom panel of Fig. [^l, 2) the range of the peak 
amplitudes is reduced, and 3) the number of free parameters is 
small. For both zones, p-mode and pseudo mode, a Lorentzian 
profile is used to model the peaks using a maximum likelihood 
estimator. It is important to note that we are interested only in the 
frequency of the centroids of the peaks, and for that a Lorentzian 
prohle is a good approximation. In Figs.|^and[^we show the re¬ 
sults of the hts for KIC 3424541 in the p-mode and pseudo mode 
regions respectively. The figures for the other stars can be found 
in Appendix A. 

After htting the two parts of the spectra, we compute the 
frequency differences of consecutive peaks: v„ - v„-i. As an ex¬ 
ample, the resultant frequency differences for KIC 3424541 are 
plotted in Fig|^ 

Two different regions are clearly visible in the frequency dif¬ 
ferences of KIC 3424541. Part of these differences - those at 
lower frequencies - are centred around 20.75 juHz correspond¬ 
ing approximately to half of the large spacing, Avnjodes/2 (see 
Table 1^, owing to the alternation between odd and even modes. 
At a certain frequency, between 1200 and 1230 pHz, the fre¬ 
quency separations increase and remain roughly constant around 
37.31 yuHz, although with a higher dispersion (see next section 
for further details). The observed acoustic cut-off frequency, Vcut, 
lies in the transition region where the differences jump from the 
p-mode region to the pseudo-mode zone. We dehne this po¬ 



Fig.3. Sections of the NAvPSD used to ht the eigenmode region 
by groups of six modes in KIC 3424541. The red line represents 
the hnal ht. 



Fig. 4. Pseudo-mode region of the NAvPSD in KIC 3424541. 
The red line is the result of the ht. 


sition as the mean frequency between the two closest points 
to this transition, represented by a red symbol in Fig. for 
KIC 3424541. The results of this analysis for the six stars studied 
here are given in Table 

3.2. Comparison between the observed and theoretical 
cut-off frequency 

As noted in Sect. [T] the cut-off frequency at the stellar surface 
scales approximately as Vcut g s/pjf^, suggesting that we 
may compare the observational cut-off frequency with that com¬ 
puted from the spectroscopic parameters. However, for our set of 
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Bmntt et al. 


Table 2. KIC number of each star, effective temperature from^ _ 

the A2Z pipeline (Mathur et al.|2010)l, Vcut, Avmodes, Avi, Av 2 , Ai, and A 2 computed as explained in the text. 


(20121 with an etTor bar of 84 K, Vmax computed using 


KIC 

Teff 

^max 

^cut 

^^modes 

Avi 

Av2 


^2 

Sun (GOLF) 

5777 

3097 + 116 

5106.41 + 61.53 

134.9 + 0.8 

70.29 + 0.65 

143.77 + 6.14 

0.018 + 0.020 

0.008 + 0.017 

Sun (VIRGO) 

5777 

3097 ± 116 

5106.41 + 61.53 

134.9 + 0.8 

71.70 + 0.56 

151.16+14.33 

0.056 + 0.038 

0.010 + 0.024 

3424541 

6180 

761.97 + 39.10 

1210.67+15.71 

38.8+1.2 

20.3 + 2.1 

39.8+1.5 

0.001+0.018 

0.011+0.016 

7799349 

5065 

568.14 + 9.87 

895.47 ± 15.38 

33.3 + 0.2 

19.3 + 2.5 

33.7 + 0.9 

0.002 + 0.017 

0.013 + 0.016 

7940546 

6240 

1099.57 ± 28.02 

1940.59 + 20.71 

58.9 + 0.4 

30.3 + 0.4 

48.9+1.0 

0.021 +0.023 

0.020 + 0.021 

9812850 

6330 

1224.38 + 52.21 

1898.08 + 24.29 

64.4 + 0.5 

31.8+ 1.6 

54.0 + 2.6 

0.009 + 0.024 

0.018 + 0.020 

11244118 

5735 

1384.07 + 50.26 

1968.69 + 21.48 

71.5 + 0.5 

37.1 + 2.0 

54.8 + 3.4 

0.014 + 0.029 

0.018 + 0.025 

11717120 

5105 

582.96 ± 5.86 

958.38+ 11.56 

37.7 + 0.3 

16.6+ 1.3 

29.9 + 0.8 

-0.003 + 0.030 

0.030 + 0.026 



Frequency (/iHz) 


Fig. 5. Consecutive frequency differences of KIC 3424541. The 
blue dashed and dotted lines are respectively the weighted mean 
value of the frequency differences in the p-mode (20.75 + 
0.20//Hz) and pseudo-mode regions (37.31 + 1.58//Hz) . The 
red symbol represents the estimated cut-off frequency, Vcut = 
1210.67+ 15.70//Hz. 


stars, the spectroscopic log g has large uncertainties (see Bruntt 
et al. 2012[) and a direct test is not possible. Alternatively, it 


has been demonstrated empirically that the frequency of max¬ 
imum mode amplitude Vmax follows a similar scale relation (see 
Bedding & Kjelds'^|2003|l. As can be seen in Fig. the ob¬ 


servational data show such a linear relation between v^ax and 
Vcut, although some stars deviate 2cr from it. In fact, from a the¬ 
oretical point of view, one does not expect this linear relation to 
be accurate enough at the level of the observational etTors. As 
shown below, values of Vcut can deviate by as much as 10% from 
the scale relation Vcut g for representative models of 

our stars. On the other hand, following [Chaplin et al. ( |2008| l, 
where theoretical values of Vmax based on stochastic mode ex¬ 
citation were computed, one can find departures from the scale 
relation as large as 25% for a model of a star with M — 1.3Mq 
and Vmax ~ 1500//Hz. 

We can gain a deeper insight by comparing the observational 
and theoretical results in the Vcm-Av plane. For the theoretical 
calculations we use a set of models computed with the CESAM 
code ( [Morel & Lebreton|2008| l. These correspond to evolution¬ 
ary tracks from the zero-age main sequence up to the red-giant 
phase, with masses between M - Q.9Mq and M = 2.QMq, and 
helium abundances between Y - 0.25 and Y = 0.28. Models 
with and without overshooting were considered. Other param¬ 
eters were fixed to standard values, such as the metallicity to 
Z = 0.02. Finally, the CEFF equation of state ([Christensen-[ 



Fig. 6. Observed cut-off frequencies Vcut versus maximum mode 
amplitude frequencies Vmax- Blue points with errors are values 
for the six stars considered in this paper and the Sun (the right¬ 
most point). A straight line through the Sun has been drawn for 
guidance. 


Dalsgaard & Dappenj 1992|l and a T(t) relation derived from a 


solar atmosphere model were used. 

In theoretical computations, a standard upper boundary con¬ 
dition for eigenmodes is to impose in the uppermost layer the 
simple adiabatic evanescent solution for an isothermal atmo¬ 
sphere. In this case eigenmodes turn to have frequencies below 
a cut-off frequency given by Wcut = 27rvcut = cjlH, where H is 
the density scale height (this is strictly speaking the cut-off fre¬ 
quency for radial oscillations, which is accurate for low-degree 
modes). To better fulfil the quasi-isothermal requirement, the 
boundary condition should be placed near the minimum temper¬ 
ature. For the Sun this corresponds to an optical depth of about 
log T = -4. Another advantage of using such a low optical depth 
is that at this position the oscillations are not too far from being 
adiabatic, at least compared to the photosphere. In our computa¬ 
tions we have used this value for all the models as the uppermost 
point, but for some red-giant stars we found that the maximum 
value of a/cut(0 in the atmosphere can be located at higher opti¬ 
cal depths; hence, we have taken that maximum value as repre¬ 
sentative of the observed Wcut- From an interpolation to the solar 
mass and radius (with Y - 0.25 fixed) we obtain from our set 
of models a value of Vcut = 5125 //Hz for the Sun, which is in 
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Fig. 7. Vcut versus Avtnodes for our set of theoretical models (grey 
dots). Blue points correspond to the observed stars, including the 
Sun (upper-right point). 



T., (K) 


Fig. 8. Relative differences 5c as a function of (black dots). 
Blue points are the same relative differences but considering the 
isothermal cut-off frequency Wa instead of Red points are 
the residuals after correcting by the polynomial fit indicated in 
the text. The fit is limited to models with Teff < 7000 K and 
Vcut > 500 yuHz. 


agreement with the observed one: Vcut = 5106.41 + 61.53pHz 
( Jimenez|2006 i. If model S ( Christensen-Dalsgaard et al.|1996| l 
is considered, the differences are a little larger, about 4%. It is 
important to remember that the solar Vcut changes with the mag¬ 
netic activity cycle ( [Jimenez et al.|20lT] l; thus, some additional 
dispersion in the observed stellar values could be due to that ef¬ 
fect. 

For the large separation, the theoretical and observed values 
can be computed in the same way. Specifically, we have com¬ 
puted the large separation Avniode.s from a polynomial fit of the 
form (v-vo)/(n-no) = 2/ where n is the radial order and 

Vo is the radial frequency closest to Vmax with radial order no- For 
the models, vo is estimated by assuming a linear relation to the 
cut-off frequency. The frequencies considered are those of the 
£ — 0 modes with radial orders n = no + 4. Pi(x) is the Legendre 
polynomial of degree i, and x is the frequency normalized to the 
interval [-1,1]. A third order polynomial has been used. From 
the Tasoul equation we expect that ao = Avmodes - Av, whereas 
other terms will mostly contain upper-layer information. We use 
only radial oscillations because, for evolved stars, the mixed 
character of some modes can introduce complications. For the 
Sun we obtain Avmodes = 134.9/tHz from the observations and 
Avmodes = 136.0pHz from model S. 

Figure shows Vcut against the large separation Avmodes for 
our set of models and the observed stars. At first glance, there 
is rough agreement between the observed and theoretical calcu¬ 
lations. but for the group of stars with Vcut ~2000 juHz it seems 
hard to explain the dispersion in their Avmodes values. 

We can proceed forward by taking into account the values 
of Teff. First, given the approximated scale relation for Vcut, we 
compute the residuals: 


. _ , _ C gVWFeff 

where v[?y( = 5079.5 yuHz is a reference cut-off frequency, corre¬ 
sponding to a model in our dataset close to the Sun (M - IMq, 
R = 1.004Ro, Y = 0.25, = 5799.2 K). The Tefr values are 

shown in Fig. As seen in the figure, the differences can be as 
large as |5c| ~ 10%. In fact, the scale relation Vcut g is 


derived theoretically by approximating the density scale height, 
H, by the pressure scale height, Hp, which is strictly valid in 
an isothermal atmosphere, and further by taking the equation of 
state of a mono-atomic ideal gas. In Fig. the blue points cor¬ 
respond to the relative differences 5a obtained by replacing in 
Eq. f^Wcut = c/2// by - cl2Hp. Thus, the first condition 
introduces the highest departures from the scale relation, which 
was expected because in the range of considered, the gas is 
mainly mono-atomic at the surface. 

Nevertheless, as seen in Fig. 5c is mainly a function of 
Teff. Hence, the scale relation can be improved if we subtract 
a polynomial fit from 5c. Only models with < 7000 K and 
Vcut > 500 pHz have been considered in that fit because this 
range includes all the stars used in the present study. For a better 
fit we also consider the dependence of 5c on Vcut- In particular, 
the red points in Fig.j^are obtained by replacing Vcut in Eq.[^by: 

>'cut = ■^'cut [l - (0.0004 - 0.0057.r - 1.0684x2 H- 0.0543y-H 

-H0.1265xy - 0.7134x2y)l (2) 

where x = - 1 and y — Vcut/v^^ - 1. The standard de¬ 

viation for these residuals is 0.3%. Hence, the modified cut-off 
frequencies, that we calibrated on the Sun deviate by Icr = 0.3% 
from the scaling relation, at least for our set of models. 

In a similar way, we compute the deviation of Avmodes from 
its expected scale relation, namely: 


^ _ (AVn,odes)° / M/Mq 

^ Avmodes V {RIRof ^ 

where (Avmode.s)'^ is the value for our reference model, the same 
as the one used in Vcut- This quantity is shown in Eig.|^for our 
set of models. 

Erom Eig. we can determine that the differences between 
the large separation computed from p-modes and the scaling 
relation can be as large as 5% for solar-like pulsators with 
Teff < 7000 K. However, as shown in this figure, the differences 
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Table 3. Values of logg derived from spectroscopy, Vmax and Vcut for each star in our sample. The mean density p is in gjcrc? and 
the large separation Avmodes in jtHz. See text for details. 


KIC 

log ftp 

loggv„„ 

lOggvcu, 

^^modes 

P 

3424541 

3.50 ± 0.08 

3.84 ±0.03 

3.86 ±0.02 

38.8 ± 1.2 

0.116 ±0.007 

7799349 

3.71 ±0.08 

3.67 ± 0.02 

3.68 ± 0.02 

33.3 ±0.2 

0.085 ± 0.001 

7940546 

4.11 ±0.08 

4.01 ± 0.02 

4.06 ± 0.02 

58.9 ±0.4 

0.270 ± 0.004 

9812850 

4.16 ±0.08 

4.06 ± 0.03 

4.06 ± 0.02 

64.4 ± 0.5 

0.327 ± 0.005 

11244118 

4.23 ± 0.08 

4.09 ± 0.03 

4.04 ± 0.02 

71.5 ±0.5 

0.388 ± 0.005 

11717120 

3.80 ±0.08 

3.69 ± 0.02 

3.71 ±0.02 

37.7 ±0.3 

0.109 ±0.002 



T„ (K) 

Fig. 9. Relative differences 6a versus T^ff for the set of models 
indicated in the text (black dots). The red points are the residuals 
after the polynomial fit. This fit is limited to models with T^ff < 
7000 K and Av> lOpHz. 



Fig. 10. Effective v*^ versus the effective for our set of 

models (grey dots) and for the observed stars (blue dots). Red 
points are observed values of the frequencies of maximum am¬ 
plitude properly scaled, (Vcut/vn,ax)GVmax- 


6A depend mainly on Teff. This was previously noted by |White1 
|et al.| ( |201 l| l and allows for a correction to Avmodes in order to 
get better a estimate of the mean density. Here we will use a 
fit similar to that considered for Vcut- As in the previous case, 
we have considered only those models with T^ff < 7000 K and 
Avmodes >10 pHz, which include all the stars used in the present 
work. The red points in Fig. correspond to the residuals ob¬ 
tained after replacing Av^ojes in Eq.j^by; 

AvLdes = [l - (0.0023 - 0.1123X- 0.7309;r2_ 

-0.0299Z - 0.0067 xz + 1.2461x^z)] (4) 

(AVmodes)*(0) 


where x = - 1 and z = Avniodes/(Avmodes)° - 1 with 

(Avmodes)*' = 135.4 juHz. The standard deviation for these resid¬ 
uals is 0.6%. Given the differences between the observed and 
theoretical values for the Sun, we have calibrated Eq. to the 
observed solar value. 

Figure 10 shows against for our set of models 

(black poinfs) and the observed stars including the Sun (blue 
points). Whereas the scattering in the theoretical values are sub¬ 
stantially reduced, the observational ones are not. Hence, if the 
error estimates are correct, we must conclude that the observa¬ 
tional cut-off frequencies do not completely agree with the the¬ 
oretical estimates. Moreover, the red points in Fig. [T^ are pro¬ 
portional to the observed values of the frequency of maximum 
amplitude scaled to the Sun, (Vcut/vmax)Gt^max- For our limited set 
of stars, Vmax follows the scaling relation better than v*^. This is a 


little surprising because, as noted before, theoretically, we would 
expect the opposite, even more so when the effective is used. 

This result can be understood in terms of the inferred surface 
gravity since we can use either v^ax or Vcut to estimate log g. The 
results are summarized in Table and compared to the spectro¬ 
scopic loggsp. To estimate the errors in loggy^,^,, we have con¬ 
sidered only the observational errors in Vmax and Teff, while for 
log^vcu, we have included the 1-cr value of 0.3% found in the 
scaling relation plus an error in derived from assuming an 
unknown composition with standard stellar values. In particular, 
we have taken ranges of T = [0.24,0.28] and Z = [0.01,0.03] for 
the helium abundance and the metallicity respectively and thus 
estimated the error in the mean molecular weight by some 
5%. As expected, the values obtained from loggy_„^,^ are quite 
similar to those reported by Bruntt et al. ( 2012| l because they are 
based on the same global parameters with very similar values. 
The values derived from v^ax and Vcut are also much closer to 
each other than those derived spectroscopically, which, as men¬ 
tioned before, have large observational uncertainties. 


4. The HIP region 

4.1. Observations 

Using Sun-as-a-star observations from the GOLF instrument on 
board SoHO, [Garcfa et al.| ( |1 998) 1 uncovered the existence of a 
sinusoidal pattern of peaks above Vcut as the result of the inter¬ 
ference between two components of a travel wave generated on 
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the front side of the Sun with a frequency v > Vcut, where the 
inward component returns to the visible side after a partial re¬ 
flection on the far side of the Sun (see Fig. 3 in Garcia et al. 
|1998| l. Therefore, the frequency spacing of ~70 //Hz found in 
the Sun corresponded to the time delay between the direct emit¬ 
ted wave and that coming from the back of the Sun (correspond¬ 
ing to waves behaving like low-degree modes), i.e. a delay of 
four times the acoustic radius of the Sun (~ 3600s). This value 
is roughly half of the large frequency spacing of the star and we 
will call it Avi. [Garcia et al.| ( |1998| l also speculated that, above 
a given frequency, a second pattern should become visible with 
a double frequency spacing (close to the large frequency sepa¬ 
ration), Av 2 ~ 140 //Hz. This pattern (see for a theoretical de¬ 
scription Kumar||1993| l - usually visible in imaged instruments 
( Duvall et al.||199f l - corresponds to the interference between 
outward emitted waves and the inward components that arrive at 
the visible side of the Sun after the refraction at the inner turning 
point (non-radial waves). 


In order to get a global estimate of the frequency of the inter¬ 
ference patterns in the HIP region of our sample of stars, two sine 
waves are htted above the cut-off frequency. The amplitudes and 
the frequencies of both sine waves are summarized in Table 
We have also re-analysed the GOLF and VIRGO data follow¬ 
ing the same procedure but at a higher frequency range, between 
7 and 8 mHz, compared to the original analyses performed by 
[Garcia et al. (19981 and Jimenez[ ( 2006[ l. In this way, we have 
been able to obtain the second periodicity at ~ 140 //Hz (see 
Tablel^. 


We tried htting various functions and opted for the simplest 
one. Indeed, the amplitude of the interference patterns decreases 
with frequency in a way close to an exponential decrease. The 
additional parameters required give more unstable hts with a 
heavy dependence on the guess parameters. Because we ob¬ 
tained the same qualitative results leading to the same classifi¬ 
cation of stars, we preferred this approach. In future studies we 
will look for a better function for the ht - i.e. one more suited to 
the observations - in a larger set of stars. We are already work¬ 
ing in this direction but this investigation is beyond of the scope 
of this paper. 


According to the fitted amplitudes, we can classify the stars 
into two groups: those where the two amplitudes of the sine 
waves, Al and Az, are similar (KIC 7940546, KIC 9812850, 
KIC 11244118), and those where A 2 is much larger than Ai 
(KIC 3424541, KIC 7799349, KIC 11717120). An example of 
each group of stars is given in Fig.[^for stars KIC 7940546 and 
KIC 11717120. It is important to notice that the error bars of the 
htted amplitudes are large because the actual amplitudes of the 
interference patterns decrease in a quasi-exponential way that 
has not been taken into account here. In this pioneering work, we 
have favoured the hts with constant amplitudes to avoid adding 
more unknown parameters to the ht. 


The stars of the hrst group have larger Avmodes than the other 
three stars, which implies a correlation with the evolutionary 
state of the stars. Starting with the Sun on the main sequence, 
the HIP pattern is dominated by interference waves partially re- 
Hected at the back of the Sun (for the solar case Ai is much 
bigger than A 2 ). When the stars evolve, the HIP pattern due to 
the interference of direct with refracted waves in the stellar inte¬ 
rior is more and more visible and becomes dominant for evolved 
RGB stars with Avmodes below ~40 //Hz. 




Fig. 11. HIPs region of KIC 11717120 (top), where the HIP pat¬ 
tern is dominated by the sine wave (red line) owing to the inter¬ 
ference between direct and refracted waves, and of KIC 7940546 
(bottom), where the two HIP patterns can be determined. The 
blue and green line in the bottom panel correspond to the two 
fitted sine waves with frequencies Avi and Av 2 respectively. The 
green and blue lines have been shifted by 0.05 for clarity. The 
red line is the actual ht. 


4.2. HIPs and stellar evolution 


In order to reproduce theoretically the periodic signals expected 
from the pseudo-mode spectrum we will follow the interpreta¬ 
tion given by [Kum^ ( 1993| l and Garcia et al. ( 1998| l and as¬ 
sume that waves are excited isotropically at a point very close to 
the photosphere. The observed spectrum of the pseudo modes is 
then interpreted as an interference pattern between outgoing and 
ingoing components, possibly including successive surface re- 
hections. In what follows we summarize the basic concepts and 
apply them to stars at different evolutionary stages. 

Let us start by using a gravito-acoustic ray theory approx¬ 
imation with the following dispersion relation (see e.g. [Gough 
|T98^|T99^ : 








—k 


(5) 


where is the radial component of the wave number, = 
L^/r^, L - l(l + 1) the horizontal component, N is the buoyancy 
frequency, is a generalized cut-off frequency given by 


ai 


2 

C 


47/2 



(6) 


and H is the density scale height. In the ray approximation the 
turning points are given by ki- - 0, whereas the ray path is 


7 
















































Jimenez, Garcia, Perez Hernandez & Mathur: Acoustic cut-oif frequency and HIPs with Kepler 


determined by the group velocity. For the spherically symmet¬ 
ric case the rays are contained in a plane with a path given by 
dOldr - vglirvy). Here, r and 0 are the usual polar coordinates 
and the group velocity Vg = (Vr, vg) - {d(jjldkr,da)ldka) has the 
following components; 


krOpc^ 2 [ to^-N^ Y 


( 7 ) 


According to the ray theory, the general solution of the wave 
equation can be expressed as a superposition of rays of the form 
t) - A^(r) exp[t'((Uf ± A;£/s)], where the two signs corre¬ 

sponds to outgoing and ingoing waves and the integral is com¬ 
puted over the ray path s, including additional reflections where 
appropriate. For every reflection, a constant phase shift must be 
introduced. In particular, for a ray travelling from an inner turn¬ 
ing point ri to an outer turning point r 2 , this integral can be ex¬ 
pressed as 


CJ 



k— dr — njA . 
Vr 


(8) 


rR 

For low-degree acoustic waves tol{2Avi^^{) ^ oj { \ jc)dr - { \ + 

{)nl2 , that is, Avoj/ ~ Av. 

Considering a wave excited very close to the surface with an 
outgoing component of amplitude and an ingoing compo¬ 
nent that emerges at the surface with an amplitude after its 
first internal traversal, and with an amplitude of A'^ [R{<jd) after a 
subsequent partial surface reflection on the back side of the star, 
the contribution to the amplitude spectrum can be expressed as: 



Fig. 12. Ray path for waves with angular degrees { - \ (red) and 
{ - 2 (green) for a 1.1 Mq model evolution sequence. For clar¬ 
ity, surface reflections have been omitted. The upper part shows 
rays with a frequency a 10% greater than the acoustic cut-off 
frequency Wa while the bottom panels corresponds to rays with a 
frequency 50% greater than Wa. The large separation indicated is 
the integral (2 l/c)“* and can be used to identified the mod¬ 

els in Fig. The leftmost panels correspond to a model near 
the ZAMS, the middle panels are for a model near the TAMS, 
and the rightmost panels are for a model in the RGB. The dotted 
circles indicate the base of the convection zone. 






'[y'. 








])■ 


( 9 ) 

where V'^f. and are visibility factors, and d is a phase 
constant that depends on the type of radiation emitted. The re¬ 
flection coefficient R is expected to be a monotonically decreas¬ 
ing function of frequency (for further details see Kumar|1993| l. 
We also assume that, for stochastically excited waves, the am¬ 
plitudes and the visibility factors will be a smooth function of 
frequency. Hence, the amplitude spectrum will be modulated 
with the periodic functions cos(w/Av^/) and cos(2(u/Av^/). 
Although in principle different values of can be expected 
in different frequency ranges and for different degrees (, their 
differences are too weak, and in the present study we fitted 
the data to just two periodic components, Av 2 = Av^j f and 
Avi = Av^//2. In addition, in some stars one of the signals could 
be masked by low-visibility factors or low reflection coefficients. 

To illustrate the problem we have considered a 1.1 Mq evolu¬ 
tion sequence and plotted in Figj^rays for two typical frequen¬ 
cies in the pseudo-mode range and three evolutionary stages. Top 
panels are for frequencies u) - l.lcoa while bottom panels are for 
a> = l.SWfl. On the other hand the star evolves from left (ZAMS) 
to right (RGB). For clarity, surface reflections have been omit¬ 
ted. 


For the model close to the ZAMS (leftmost panels) the in¬ 
ward rays emerge at the surface making angles between 160° 
and 180° with the outward component. In this case, the outward 
and inward components of a given wave can be simultaneously 
visible only if they lie too close to the limb. Hence, their sig¬ 
nal corresponding to Av 2 - Av will be highly attenuated in the 
power spectrum. On the other hand, after one surface reflection 
the angle between the outward and inward components lies in 


the range 0° to 35°. This interference would give the Avi ^ Av/2 
separation and, hence, had it the same intrinsic amplitude as the 
former, the Avi separation would be easier to observe. But in 
this case the inward rays are only partially reflected on the far 
side (e.g. |Kunw ( |1993| l) and hence the attenuation factor, corre¬ 
sponding to R(cij) in Eq. [9j is high. Since the actual observations 
show that in the Sun the dominant separation in the HIP pattern 
is Avi ( |Garcfa et al.|199^ we can use this case as a qualitative 
reference between the two competing factors. For model S, we 
obtain Avi =72 /iHz with our fit, for which a frequency range 
~ [5200,6600] pHz was used. This result is in agreement with 


the observational value of 70.46 + 2 jAiz found by Garcia et al. 
( 1998| l. As mentioned in section 4.1, we re-analysed the solar 
data (GOLF and VIRGO/SPM) and obtained Avi~ 70 pHz (see 
Table ig, in agreement with previous values given in the litera¬ 
ture. At higher frequencies we re-fitted the data with two sinu¬ 
soidal components and found Av2~140jL(Hz (Table|^. This value 
is in agreement with the theoretical prediction already described 
by Garcia et al. (1998) if we assume that the photosphere is the 
only source of partial wave reflection. However, in this case the 
amplitude is much smaller because only waves close to the limb 
contribute to this interference pattern. 

As the stars evolve, the angle between the outward ray and 
the first surface appearance of the inward refracted ray becomes 
smaller for the { - I pseudo-modes. This is clearly apparent in 
Fig. 12 Hence, the Av 2 = Av interference pattern becomes easier 


to observe (higher observed A 2 amplitudes). This phenomenon 
is due to the transition between acoustic rays smoothly bending 
through the stellar interior, and waves where the density scale 
height in the stellar inner structure becomes close to their wave¬ 
length and the ray behaviour at the inner turning point is closer 
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to a two-layer reflection. Finally, in the most evolved model in 
Fig. (upper and bottom right panels) the { - \ waves have 
a gravity character close to the centre that bends the ray into a 
loop. 

In order to reproduce the observed power spectrum schemat¬ 
ically we used Eq. and considered ray traces for a contin¬ 
uum spectrum of waves with frequencies between the acous¬ 
tic cut-off frequency at the surface and about a 1.5 times that 
value. This frequency interval approximately spans the observed 
pseudo-mode frequency range. Since we are interesting only in 
reproducing the periodicities of the pseudo-mode spectra, con¬ 
stant amplitudes were considered, and visibility factors either 
constant or proportional to cos 0/2, where 0, is the angular dis¬ 
tance between the outward and inward waves. We determine the 
periodic signal in the spectrum by a non-linear fit to the above 
equation. 

Fig. 13 shows the angular distance between the inward and 
outward components before any surface reflection against the 
large separation Av for models in a 1.1 Mq, Y - 0.28 evolu¬ 


tion sequence with the code and physics indicated in section 3.2 


Here, mean values of the angular distance are computed by av¬ 
eraging those of the ray traces in the frequency range indicated 
above. Pseudo-modes with ( - \ and 2 are averaged separately. 
Since no surface reflection is considered, these waves form the 
Av 2 pattern. 

Looking at the behaviour of the ( - \ pseudo-modes in 
Fig|T^(red points) we may conclude that for M - 1.1 Mq stars 
that evolve to a Av below some given threshold between 90 //Hz 
and 60 //Hz (corresponding to evolved main sequence stars as the 
one shown in the middle panel of Fig.[T2|), the Av 2 pattern should 
become visible. Indeed, a smaller angular distance means more 
disc-centred interferences favouring a higher A 2 . The critical Av 
slightly changes with mass, being lower for higher masses. In 
any case all the stars in Table except the Sun, are evolved to a 
point where the pattern with Av 2 could be easily observed, as is 
in fact the case. 

Regarding the { - 2 pseudo-modes, it can be seen in Fig. 
that the angle distance always remains above some 150° and 
hence this degree contributes only to the Avi period in the 
pseudo-mode spectrum. The higher dispersion in the angular dis¬ 
tance for models with a large separation between Av = 60 and 
90//Hz (evolved main sequence stars and subgigants), clearly 
visible in Fig. deserves a comment. According to the dis¬ 
persion relation used in this study, for models evolved near the 
TAMS the characteristic frequencies of the ( - 2 waves become 
complex in a small radius interval above the core, thus increas¬ 
ing the resonant cavity of the pseudo-modes within a limited 
frequency range. An example of the ray path of these kinds of 
waves is shown at the bottom middle panel of Fig. (green 
line). Here, the ray trace close to the centre is not of the acoustic 
type and of the emerging ray is consequently scattered. Although 
these results may be questionable in terms of the validity of the 
approximation used here, they do not have any observable con¬ 
sequences. 

For the three most evolved stars in our dataset the interfer¬ 
ence pattern coming from waves reflected back at the surface 
with period Avi is not observed. Because, as noted before, for 
evolved stars the signal from the ( - \ pseudo-modes do not 
suffer any attenuation from surface reflections, it is possible that 
this explains why for stars evolved to the RGB the Avi signal 
becomes completely masked. However, it is interesting to note 
the another, possible superimposed, cause for this circumstance. 
In principle, one might expect that radial waves propagates all 
the way from the surface to the centre and hence contributes to 



Av( pHz) 


Fig. 13. Angular distance between inward and outward waves 
against Av for a 1.1 Mq evolution sequence. Frequencies be¬ 
tween the acoustic cut-off frequency at the surface and about 
50% greater than o/a are considered. For every model in the se¬ 
quence, frequency-averaged values of the angular distance and 
Icr dispersions are calculated for € - \ (red points) and t - 2 
(green points) waves separately. An angle of 180° means no in¬ 
ner reflection. 


the signal only once reflected back, thence with a period Avi. 
However there is a point in the evolution where the cut-off fre¬ 
quency in the core rises above the typical pseudo-mode frequen¬ 
cies, in which case the ray theory introduces an inner reflec¬ 
tion. Although a piano-parallel approximation is questionable in 
terms of the wavelength-radius relation when we are too close 
to the centre, this approximation allows us to estimate the trans¬ 
mission coefficient T, as for the one dimensional p roblem[^If the 
equation T - 1 /(I - 1 - e^^) is used with K - 2 
ri being the point where kr - 0, it happens that, for red gi¬ 
ants stars, radial waves in the observed frequency range of the 
pseudo-mode spectrum are mostly reflected at the core edge. 
Thus at this stage radial pseudo-modes also contribute to the Av 2 
signal and hence the Avi period would hardly be observed. 

Let us now discuss the relation between the large separation 
obtained from the eigenmodes, Avmodes and those corresponding 
to the pseudo-modes, Avi and Av 2 . First, we have verified that 
for waves with { - 0 the phase travel time is always close to 


the asymptotic acoustic value, f kdl ^ I jc dr. Hence, when 
these acoustic waves contribute to the HIPs after a surface reflec¬ 
tion they give Avi ^ Avmodes/2. In addition, most of the £ - 2 
and, in some stars, € - \ travelling waves have similar acoustic 
characteristics, which also contributed to the same interference 
period. This is in agreement with Table where both frequency 
separations match within the errors. 

The case of Av 2 is different. The grey points in Fig. 14 are the 
periods of the interference pattern computed according to Eq. 
for travelling waves and no surface reflection. A wide range of 
masses and evolutionary stages from main sequence to the base 


' For radial oscillations the full adiabatic equations are of second or¬ 
der and hence a WKB analysis can be done without the approximations 
assumed in the dispersion relation Eq. jg. The qualitative results given 
in this paragraph rely solely on the assumption of a one-dimensional 
problem. 
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of the red-giant branch were included. The red points are the ob¬ 
served Av 2 values while blue points are 2Avi when observed, 
including the Sun. The black line corresponds to Avmodes = Av 2 . 
From Fig. we can see that for Avmodes = 60-80/tHz there are 
models with Av 2 < Avmodes- In fact they correspond to evolved 
stars up to the end of the main sequence or the subgiant phase, as 
it is the case for our sample of stars. Although the order of mag¬ 
nitude of the phase delay found for our models is similar to the 
observed one for this type of star, it is also apparent from the hg- 
ure that KIC 1124411 (Avmodes ~ 70 pWz) has a large separation 
- too high for this delay to appear. We obtain a mass of ~ IMq 
either by applying the scaling relations to Vmax and Avniode.s by 
doing an isochrone ht. The upper right corner of Fig|^is repre¬ 
sentative of this star. Hence we expect to observe the signal Av 2 
from waves reflected in the interior but our computations do not 
shown any signihcant delay in the phase of the € - \ compared 
to the radial oscillations. Further work is need here; in particular, 
the asymptotic theory that we have used could be inadequate for 
such details. 

For the stars with the lower Avjnodes two of them have Av 2 - 
Avmodes and one star (KIC 11717120) has Av 2 < Avmodes- With 
our isochrones htting, we hnd that the latter is at the base of the 
RGB as well as KIC 7799349, thus we have two different results 
for stars with very similar parameters. With our simple simula¬ 
tion for models at this evolutionary stage we have both signals 
since, as noted above, radial oscillations are partially reflected at 
the core edge and for these stars, Av 2 - Avmodes- The observed 
period will be a weighted average, but for a better comparison 
with the observations proper amplitudes need to be computed. 



Fig. 14. Grey points are periods of the pseudo-mode interference 
pattern for evolution sequences with masses between Q.9Mq and 
2Mq. For the sake of clarity only waves with no surface re¬ 
flections are shown. Red points correspond to the observed Av 2 
values and blue points correspond to 2Avi when observed. The 
black line corresponds to Avmodes = Av 2 


5. Conclusions 

In this study the acoustic cut-off frequency and the characteris¬ 
tics of the HIPs have been measured in six star with solar-like 


pulsations observed by the Kepler mission. A comparison with 
the observed Vmax shows a linear trend with all stars lying in 
~2-cr. As a result, the values of log g derived from Vcut agree 
to within the same accuracy with those derived from Vmax but 
are substantially different in some cases from the values derived 
from spectroscopic hts (see Table |^. When comparing the stars 
with the models in the Vcut-Av plane we found a departure in Vcut 
from the expected values. While using theoretical information it 
is possible to calculate a measurable function v*^ = Vcut) 

that scales as g/ with no more than 0.3% deviation (for 
our set of models representative of our star sample and the Sun), 
we hnd that the observational v*^ does not follow this relation 
so accurately. Rather surprisingly, we found that the frequency 
of maximum power, Vnax. follows the linear relation within er¬ 
rors. Hence, for the evolved stars considered in the present study 
it must be concluded that gives an even better estimate of 
logg than Vcut with the current uncertainties. 

Given the observed characteristics of the HIPs, our set 
of six stars can be divided into two groups. In three stars: 
KIC 3424541, KIC 7799349, and KIC 11717120, the only vis¬ 
ible pattern is the one due to the interference with inner re¬ 
fracted waves in the visible disc of the star and not close to the 
limb. In the other three stars ( KIC 7940546, KIC 9812850, and 
KIC 11244118) the pattern due to partial reflection of the inward 
waves at the back of the star is also detected. This different be¬ 
haviour is related to the wide spacing of the star, which reveals 
a dependence on the evolutionary stage. 

When present, the period Avi, which corresponds to the in¬ 
terference of outward waves with their inward counterpart once 
reflected on the far side, always agrees with half the large sep¬ 
aration, although, like the Sun, their values are a little higher 
(this result is also found theoretically, the large separation de¬ 
rived from the pseudo-modes being closer to f dric). However, 
the period Av 2 corresponding to waves refracted in the inner part 
of the star can be substantially different from the large frequency 
separation. We interpret this by claiming the presence of waves 
of a mixed nature. For these stars, the phase f kdl is smaller 
compared to the radial acoustic case. Although a detailed anal¬ 
ysis is beyond the scope of the present study, simple ray theory 
calculations reveal that such a phenomenon is expected with the 
same order of magnitude. 

The pseudo-mode spectrum reveals information not only 
from the surface properties of the stars through Vcut but also from 
the interior, at least in the cases where the period Av 2 is lower 
than the large separation derived from the eigenmodes. We are 
aware that the interference periods derived observationally are 
power-weighted averages and for a proper comparison with the¬ 
oretical expectation some work along these lines should be ad¬ 
dressed. This might be accomplished by computing transmission 
coefficients so that more realistic theoretical simulations of the 
interference phenomenon can be done. 
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Appendix A: Figures of the other stars 



Fig. A.l. Fitted (red line) NAvPSD of KIC 7799349 by groups 
of nine peaks in the low-frequency range (p modes). 



Fig. A.2. Fitted (red line) NAvPSD of KIC 7799349 in the high- 
frequency range (pseudo modes). 



Fig.A.3. Consecutive frequency differences (KIC 7799349); 
that is, the separations between the fitted peaks for NAvPSD. 
Two levels appear, one around Av/2 corresponding to p-modes 
(blue dashed line with a weighted mean of 16.64 + 0.09 yuHz) 
and other around Av corresponding to pseudo-modes (blue dot¬ 
ted line with a weighted mean of 35.23 + 1.52 pHz). The red 
symbol is the estimate of the acoustic cut-off frequency (895.47 
+ 15.38 pMz). 



Fig. A.4. Fitted (red line) NavPSD of KIC 7940546 by groups of 
five peaks in the low-frequency range (p modes). 
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Fig. A.5. Fitted (red line) NAvPSD of KIC 7940546 in the high- 
frequency range (pseudo-modes). 
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dashed line with a weighted mean of 29.50 ± 0.07 //Hz) and 
other around Av corresponding to pseudo-modes (blue dotted 
line with a weighted mean of 45.71 ± 2.12 //Hz). The red sym¬ 
bol is the estimate of the acoustic cut-off frequency (1940.59+ 
20.71 //Hz). 
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Jimenez, Garcia, Perez Hernandez & Mathur: Acoustic cut-off frequency and HIPs with Kepler 
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Fig.A.9. Consecutive frequency differences (KIC 9812850); 
that is, the separations between the fitted peaks for NAvPSD. 
Two levels appear, one around Av/2 corresponding to p-modes 
(blue dashed line with a weighted mean of 32.53 + 0.12 //Hz) 
and other around Av corresponding to pseudo-modes (blue dot¬ 
ted line with a weighted mean of 58.13 + 2.39 //Hz). The red 
symbol is the estimate of the acoustic cut-off frequency (1898.08 
+ 24.29 //Hz). 



Fig.A.ll. Fitted (red line) NAvPSD of KIC 11244118 in the 
high-frequency range (pseudo modes) 
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Fig.A.12. Consecutive frequency differences (KIC 11244118); 
that is, the separations between the htted peaks for NAvPSD. 
Two levels appear, one around Av/2 corresponding to p-modes 
(blue dashed line with a weighted mean of 35.54 + 0.08 //Hz) 
and other around Av corresponding to pseudo-modes (blue dot¬ 
ted line with a weighted mean of 61.02 + 2.19 //Hz). The 
red symbol is the estimation of the acoustic cut-off frequency 
(1968.69 +21.48 //Hz). 


Fig. A.IO. Fitted (red line) NAvPSD of KIC 11244118 by groups 
of five peaks in the low-frequency range (p modes). 
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Jimenez, Garcia, Perez Hernandez & Mathur: Acoustic cut-oif frequency and HIPs with Kepler 
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Fig. A.13. Fitted (red line) NAvPSD of KIC 11717120 by groups 
of sixteen peaks in the low-frequency range (p-modes). 



Fig.A.15. Consecutive frequency differences (KIC 11717120); 
that is, the separations between the htted peaks for NAvPSD. 
Two levels appear, one around Av/2 corresponding to p-modes 
(blue dashed line with a weighted mean of 16.71 + 0.09 pHz) 
and other around Av corresponding to pseudo-modes (blue dot¬ 
ted line with a weighted mean of 29.85 + 0.66 yuHz). The 
red symbol is the estimation of the acoustic cut-off frequency 
(958.38 + 11.56 a(Hz). 



Fig.A.14. Fitted (red line) NAvPSD of KIC 11717120 in the 
high-frequency range (pseudo-modes). 


14 






























